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Abstract 

We analyze the tensor mode perturbations of static, spherically symmetric solutions of the 
Einstein equations with a quadratic Gauss-Bonnet term in dimension D > 4. We show that 
the evolution equations for this type of perturbations can be cast in a Regge- Wheeler- Zerilli form, 
and obtain the exact potential for the corresponding Schrodinger-like stability equation. As an 
immediate application we prove that for fl / 6 and a > 0, the sign choice for the Gauss-Bonnet 
coefficient suggested by string theory, all positive mass black holes of this type are stable. In the 
exceptional case D = 6, we find a range of parameters where positive mass asymptotically flat 
black holes, with regular horizon, are unstable. This feature is found also in general for a < 0. 
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I. INTRODUCTION 



Alternative gravity theories in higher dimensions have been attracting considerable at- 
tention, particularly the Einstein-Gauss-Bonnet (EGB) theory, which emerges as the low 
energy limit of string theory. The EGB lagrangian is a linear combination of Euler densities 
continued from lower dimensions. It gives equations involving up to second order derivatives 
of the metric, and has the same degrees of freedom as ordinary Einstein theory. A particular 
choice of the coefficients in front of the Euler densities gives theories where the local Lorentz 
symmetry is enlarged to a local (A)dS symmetry [1, 2]. Interesting solutions to the EGB 
equations, many of them relevant to the development of the AdS — CFT correspondence 
[3], include a variety of black holes in asymptotically Euclidean or (A)dS spacetimes. These 
solutions could be found mostly because they are highly symmetric. Analyzing their lin- 
ear stability, however, confronts us with the high complexity of the EGB equations, since 
the perturbative terms break the simplifying symmetries of the background metric. In this 
Letter we report on the stability of spherically symmetric, static solutions of the quadratic 
EGB theory. These are preliminary results of ongoing work on the stability of EGB black 
holes with arbitrary Einstein manifolds as horizons [4] 



II. TENSOR PERTURBATIONS OF SPHERICALLY SYMMETRIC EGB SPACE- 
TIMES 

The lowest order Einstein-Gauss-Bonnet (EGB) vacuum equations are 

= Q b a = AG (0)b a + G {1)b a + aG {2)b a (1) 

Here A is the cosmological constant, G^) ab = g a b the spacetime metric, G(\) ab = R a b — \Rg a b 
the Einstein tensor and 

f~i a p de p ca r) p c p da opcpa,ppa ^ ro / p ef p cd a p d p c i p2\ 

"(2)6 — -n-cb t^de ~ ^tid ticb ~ ^tib tic + tXttb — —0 b [K c d tief ~ 4r£ c tid + ti J 

(2) 

the quadratic Gauss-Bonnet tensor. These are the first in a tower G^ b a , s = 0, 1, 2, 3, ... of 
tensors of order s in R a b cd constructed by Lovelock [5]. As shown in [5], the most general 
rank two, divergence free symmetric tensor that can be constructed out of the metric and 
its first two derivatives in a spacetime of dimension d, is a linear combination of the G( s ) b a 
with 2s <d [5]. 

Here we consider the spherically symmetric case, a spacetime of dimension D = n + 2 with 
metric 

ds 2 = -f(r)dt 2 + jT-rdr 2 + r^g^dx^ (3) 

where gijdx l dx^ is the line element of the unit n— dimensional sphere S n . We use indices 
i, j, k, I, to, ... and a bar for tensors and operators on S n , whereas a,b,c,d,... are generic 
indices. The nonzero components of the Riemann tensor of the metric (3) are 

Rtr r = - y, Rit jt = Rir jT = V = {^f) - ^ ( 4 ) 
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Inserting (4) in (1) we find that (3) solves the EGB equation if 



/(r) = l-r 2 V(r), (5) 

and ip(r) satisfies [6] 

an(n — l)(n — 2) , . . 9 n , . . A a 

— ^ ^ -i){rf + -Mr) = 6 

4 ^ K J 2^ K J n+1 r n+l V ; 

We consider tensor perturbations around (3) 

9ab -> 9ab + Kb- (7) 

which are those satisfying h ab = unless (a, b) = Tensor perturbations are believed to 

be the only potentially unstable modes in ordinary Einstein theory [7]. We choose the gauge 
where h a b is transverse traceless. This is easily seen to imply that the restriction of h a b to 
the sphere is transverse traceless, and so can be expanded using a basis of eigentensors of 
the laplacian [8]. Thus, we need only consider the case 

hij(t,r,x) = r 2 (f)(r,t)hij(x) (8) 

where 

V'T/,//,, ~h, r V7/, ; 0. //"7/, ; (9) 

Solutions to equations (9) are worked out in [8], where it is shown that the spectrum of 
eigenvalues is 7 = — l(l + n — 1) + 2, I = 2,3, 4, ... The components of the first order variations 
5G( s ) b a , s = 0,1,2 under (8) are trivial unless (a, b) = After a long calculation the 

components are found to be 

5G {0)i j = (10) 



SG {1)i j = 5R { 

and 



w (11) 



*W = { (0 - f<t>") (j^j) l-rf + (n - 3)(1 - /)] 

+0' {^f) [(n - 3) ((n - 2)(/ 2 - /) - rf) + r 2 (f 2 + /"/) + (3n - 7)r/'/] 

+0 ("a^) [ r2/ " + 2(n " 3)r/ ' + (n " 3)(n " 4)(/ " 1} ] } ^ (12) 

Setting <f>(r,t) = e ut x(r) the linearized EGB equations 

5G (1)a b + «5G (2)a fe = (13) 
around the metric (3) reduce to a second order ODE for x( r )- By further introducing, 

$(r) = x(r)if(r) (14) 
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with, 



K(r) = W 2 " 1 ^ 2 + a{n - 2) Jfa - 3)(1 - f) - r^j (15) 

and switching to "tortoise" coordinate r*, defined by dr*/dr = 1//, this ODE can be cast 
in the Schrodinger form, 

' '" + V(r(r*))$ = -co 2 ® = (16) 



dr* 

The solutions will therefore be stable if (16) has no negative eigenvalues. On the other 
hand, a negative eigenvalue (E < 0) signals the possibility of an instability that requires 
also consideration of the normalization of the corresponding eigenfunctions (see, e.g. [7] for 
details). 

The explicit form of the potential V(r) as a function of r and the parameters of the theory 
is rather lengthy. We notice however that if we introduce the function, 

/ /(2- 7 ) \ / (1 - af»y + q(n - 3) [(n - 4) (1 - /) - 2rf>] \ 
q \ J\ r2 + a ( n _ 2 )[(n-3)(l-/)-r/1 J 1 J 

the potential is given by, 

V(r) = q(r) + (f± In (*)) ' + / ^ (/ ^ In (K)) (18) 

Eqns (14)-(18) are the main result of this paper, (18) being the exact potential of the 
Schrodinger-like stability equation for spherically symmetric EGB blackholes of arbitrary 
mass and cosmological constant. Clearly, it can be applied to the cosmological solutions 
of the EGB equations that result by setting // = 0. Moreover, our results are readily seen 
to reproduce those in [7] in the a = (Einstein gravity) limit, which was also extensively 
studied by Kodama and Kodama and Ishibashi (see, e.g., [9] and references therein), as well 
as the restricted cases studied in [10] and [11]. In what follows, as an application of the 
formalism, we analyze briefly the case A = for general n, and also the n = 3 and n = 4 
BTZ black holes [2] . The general case will be considered in a more extended version of this 
paper, currently in preparation [4]. 



III. STABILITY OF EINSTEIN-GAUSS-BONNET BLACK HOLES 

We recall that for A = 0, on account of (5) and (6), for asymptotically flat Einstein- 
Gauss-Bonnet black holes with regular horizon /(r) takes the form [6], 



f(r) = 1 + 



n 



l)(n-2)a 



1 + 



4(n - l)(n - 2)a/x 



nr 



n+l 



(19) 



where \i > corresponds to positive mass. We consider first a > which is the relevant 
case for string motivated theories. Then, for any \x > 0, there is a regular horizon at r = 
and /(r) grows monotonically from zero to one as r grows from r H to infinity. From (19), 



fx = n[a(n — l)(n — 2) + 2r H ]r^ 



(n-3) 



(20) 
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Going back to (16), a sufficient criterion for stability is that V(r) is positive for r > rn- 
If we consider (18), we notice that the second term on the R.H.S. is positive definite in 
all cases, while, a long computation shows that the first and third terms are also positive 
definite for r > r H , for n = 3 and all n > 4, so all these cases are stable under tensor 
perturbations. The n = 4 case is exceptional. Here we notice that, since V(r(r*)) is 
bounded in — oo < r* < +oo, with V(r(r*)) — > for r* — > ±oo, a sufficient condition for 
the existence of a bound state of negative energy is [12], 



j — < 



V(r(r*))dr* < 
This can be written as an integral over r, 



(21) 



r+oo 

y < 

Jr H 



{V{r)/f{r))dr < 



(22) 



The second term on the right in (18), divided by /, is positive, while the third, divided by 
f(r), is a total derivative in r, and gives a vanishing contribution on account of its behaviour 
for r — > r# and r — > +oo, as is easily seen from (15). The "dangerous" contribution comes 
then from q(r)/f(r). In fact, since q(r) contains the (positive) factor (2 — 7), which can be 
arbitrarily large for spherical horizons, while the other terms in (18) are independent of 7, 
the condition, 



+00 



g(r)/((2- 7 )/(r))dr<0 



(23) 
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implies that (21) will be satisfied for a sufficiently large 7. Note that a has dimension r 2 
and that for n — 4, /i has dimension r 3 . Introducing z = fia~ 3 ^ 2 in (20) we find 



(2z + 2 v / 16 + z 2 ) 2/3 -4 
(2z + 2V16 + 2 2 ) 



1/3 



(24) 



so that rn 
gives 



+ as z — > + (/x — > + ). Setting n = 4 in (17) and defining x = r(/ia) 



-1/5 



(2-7)/ 



(no) 



-2/5 



2(x 5 + 6) 2 - 75 
2x 2 (x 5 + 6)(x 5 + 1) 



(25) 



The integral in (23) can be given in closed form using (25), but the expression is too long 
and difficult to handle. We may however show that the integral in (23) is negative if we first 
change variables to u — 1/x 



r 



q{r) 



dr = 



(2- 7 )/(r) (/ia)V5 



,10 



2(1 + 6m 5 ) 2 - 75m 
2(1 + 6w 5 )(l + u 5 ) 



du, 



(26) 



and then note from (24) that Xh — /i 4//5 o; _6 / 5 /12 for /i > 0, which implies that the upper 
limit of the R.H.S. integral above tends to infinity as /i — > + . Since the integrand stabilizes 
in —1/4 for large u, the integral is certainly negative for sufficiently small fi. To illustrate 
this point we display in Figure 1 the potential V(r(r*)) as a function of r* for two a = 1 
cases. Figure l.a shows the potential corresponding to a small \i case where (21) holds, 
whereas Figure l.b shows the potential of a large ji solution. This is positive definite and 
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therefore does not allow bound states. 

In closing this section we remark that, in space time dimensions D = 5 and D = 6, EGB 
black holes with a cosmological constant contain as particular cases the corresponding BTZ 
black holes [2]. In the notation of this Letter and that of [2] we have, 

a = f/2 , A = , n = 3£ 2 {M + l)/4 , (for D = 5) 

a = f/Q , A = -5/£ 2 , n = 2£ 2 M , (for D = Q) (27) 

Interestingly, we find that all D = 5 solutions are stable, while all solutions are unstable 
for D = 6. We recall that these cases were actually excluded in the analysis in [2], on 
considerations based on cosmic censorship. 

IV. COMMENTS AND CONCLUSIONS 

Summarizing the results reported in this Letter, we have found an explicit form for the 
Schrddinger-like equation governing the evolution of linear tensor perturbations of static 
spherically symmetric solutions of EGB vacuum equations. As a first application we proved 
the stability of (asymptotically flat) EGB black holes with positive mass and coupling con- 
stant a, in dimension D = n + 2, for n = 3, and n > 4. In the case n = 4 we found the 
unexpected result that the EGB black holes are stable only for sufficiently large mass. The 
nature of the instability of the small mass black holes is an intriguing question, outside the 
scope of the present work (a thermodynamic instability of some asymptotically (A)dS EGB 
black holes was also found in [13]). Preliminary results indicate that in the a < case, for 
all n > 3 there are solutions that represent static black holes with regular horizons, that 
are, however, unstable under tensor perturbations. The results obtained in this Letter are 
straightforwardly extended to blackholes with non positive constant curvature horizons, as 
those studied in [14]. These are currently being analyzed together with other black holes 
having more general manifolds as horizons. 
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FIG. 1: The potential V(r(r*)) as a function of r*. Figure 1 a) corresponds to fi = 0.4, a = 1, 
while for Figure 1 b) we have taken /i = 8, a = 1. The values of the integral in (21) are —9.008... 
for Figure 1 a), and +10.39... for Figure 1 b). 
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